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^ ] Abstract 

We study a two-spin quantum Turing architecture, in which discrete local rotations 
£NJ | {a m } of the Turing head spin alternate with quantum controlled NOT-operations. 

Substitution sequences are known to underlie aperiodic structures. We show that 
| ■ parameter inputs {a m } described by such sequences can lead here to a quantum 

dynamics, intermediate between the regular and the chaotic variant. Exponential 
| parameter sensitivity characterizing chaotic quantum Turing machines turns out to 

be an adequate criterion for induced quantum chaos in a quantum network. 

> : 

on ■ 1 Introduction 

Models described by one-dimensional Schrodinger equations with quasi-periodic potentials 
display interesting spectra: this kind of potential is intermediate between periodic ones, 
leading to energy bands and extended states, and truly random potentials, which cause 
localisation A super-lattice e.g. may be made of two species of doped semiconductors 
producing a one-dimensional chain of quantum wells. A qualitative model to describe the 
corresponding wave-functions is given by the tight-binding model, Hip(n) = if)(n+l)+if>(n— 
1) + \V(n)il)(n) , if) £ / 2 (Z) , where V(n) represents the effect of quantum well n, and A is 
a positive parameter playing the role of a coupling constant 0. An interesting description 
£> . results when this super-lattice is constructed by means of a deterministic rule. The simplest 
rule obtains when the two species alternate in a periodic way. But in general the rule 
will be aperiodic. One widely studied example is the Fibonacci sequence, which is quasi- 
periodic ||: given two 'letters' a and b, one substitutes a — > £(a) = ab and b — > £(&) = a. 
Iterating this rule on a one thus generates the sequence abaababaabaababa • • • , in which 
the frequency of a's is given by the golden mean (5 1 / 2 — l)/2. Other examples of such 
substitution rules are the Thue-Morse sequence (non quasi-periodic), which is obtained 
through the substitution a — > £(a) = ab , b — > £(&) = ba giving abbabaabbaababba • • • , 
and the period-doubling sequence (non quasi-periodic), a — >• £(a) = ab , b — > £(&) = aa 
P|. Accordingly, a piece-wise constant potential {V n ;n £ Z} based on such such a rule, 
e.g. V\ — a, V2 — b, V3 = a, ■ • ■ , is called a 'substitution potential'. The potential of the 
Fibonacci sequence has one-dimensional quasi- crystalline properties. 

In recent years problems of quantum computing (QC) and information processing have 
received increasing attention. To solve certain classes of problems in a potentially very 
powerful way, one tries to utilize in QC the quantum-mechanical superposition principle 
and the (non-classical) entanglement. In the models of QC based on quantum Turing 
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machines (QTM) |3|, [|, the computation is characterized by sequences of unitary trans- 
formations (i.e. by the corresponding Hamiltonians H acting during finite time interval 
steps). Benioff has studied the tight-binding model in a generalized QTM ||, where QC 
is associated with different potentials at different steps as an 'environmental effect'. These 
kinds of influences may introduce deterministic disorder which would degrade performance 
by causing reflections at various steps and decay of the transmitted component ||. 

Here we investigate an iterative map on qubits which can be interpreted as a QTM 
architecture J7|: local transformations of the Turing head controlled by a sequence of ro- 
tation angles {a m } , m = 1, 2, • • • (parameter inputs) alternate with a quantum-controlled 
NOT-operation (QCNOT) with a second spin on the Turing tape. Those angles at steps 
n = 2m — 1 are reminiscent of the potentials V m introduced before. In the present paper 
we will investigate the Fibonacci rule and the Thue-Morse case mainly with respect to the 
local dynamics of the Turing head, which will be shown to reflect the degree of 'random- 
ness' of the substitution sequences. The various types of aperiodic structures have been 
characterized up to now by the nature of their Fourier spectra only 0. 

2 Quantum Turing machine driven by substitution se- 
quences 

The quantum network || to be considered here is composed of two pseudo-spins \pf^; p = 
— 1, 1; // = S, 1 (Turing-head S, Turing-tape spin 1, see figure H) so that its network state 
l^) lives in the four- dimensional Hilbert space spanned by the product wave-functions 
\j^k^) = \jk). Correspondingly, any (unitary) network operator can be expanded as a 
sum of product operators. The latter may be based on the Si7(2)-generators, the Pauli 
matrices (jj, j = 1, 2, 3, together with the unit operator 1^. 

The initial state \ipo) will be taken to be a product of the Turing-head and tape wave- 
functions. For the discretized dynamical description of the QTM we identify the unitary 
operators U n , n = 1, 2, 3, • • • (step number) with the local unitary transformation on the 
Turing head S, Ui s J, and the QCNOT on (S, 1), U^' 1 ^, respectively, as follows: 

U 2m -i = exp (-ia[ S) a m /2 s j (1) 

u 2m = = Pi5_ l + A ( ? = {u {S > l) y , (2) 

where = \jf s ^ s \j\ is a (local) projection operator, and the Turing head is externally 
driven by substitution sequences {a m } specified by ai,a 2 - Here we restrict ourselves to 
the quasi-periodic Fibonacci - (qf) and Thue-Morse sequence (tm), respectively: 

qf qf qf qf qf , s 

Oil = Oil, a 2 = a 2 , «3 — Oil j «4 — Oil , Olg = «2 , • • • (o) 

a\ m = ai , al m = a 2 , a* 111 = a 2 , a* 111 = an , a^ 111 = a 2 , ■ • • . (4) 

First, we consider the reduced state-space dynamics of the head S and tape-spin 1, 
respectively, 



af\n) = Tr(p(f)af) = ty n \af ) ® i«|V> n > 
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a£\n) = Tr = ® &$>\i/> n ) , (5) 

where \ip n ) is the total network state at step n, and Oj (n) are the respective Bloch- vectors. 
Due to the entanglement between the head and tape, both will, in general, appear to be 
in a 'mixed-state', which means that the length of the Bloch- vectors in (||) is less than 
1. However, for specific initial states the state of head and tape will remain pure: 
As \±p = ^ ± are the eigenstates of a{ 1] with &P\±p = ±|±) (1) , the 

QCNOT-operation [A 5 ' 1 ) of equation (^) cannot create any entanglement, irrespective of 
the head state \(ff s \ i.e. 

fj&l) g, | +) (D = | +) (1) 

^l/^l-f = (6) 

As a consequence, the state \ip n ) remains a product state for any step n and initial product 
state \if)f) = \ipo) (s) (g) |±) (1) with |y? ) (5) = exp (-zafVo/2) |- 1) (5) , so that the Turing 
head performs a pure-state trajectory ('primitive', see J7|) on the Bloch-circle (a[ S \n) = 

i^> = \^) {S) ® i±) (1) , (4 5) n±)) 2 + w±)) 2 = i . 

Here erj '(n|±) denotes the Bloch- vector of the Turing head S conditioned by the initial 
state \tpo)- From the Fibonacci sequence ([|) and the property ([]) it is found for \(f^)^ <S> 
|+)W,n = 2m, and = ct = that 

4 5) (2m|+) = sinC 2m (+) , af ) (2m|+) = -cosC 2m (+) , (7) 

and ct^, (2m — 1|+) = erf, (2m |+) , where the cumulative rotation angle is 

m r 

C 2m (+) = ^2 a j = aim + (a 2 -a 1 )m / , 

3=1 

with w! (< m) being the total number of angles ct 2 up to step 2m. For the cumulative 
rotation angle C n (—) up to step n we utilize the following recursion relations 

C 2m (— ) = — C 2m _i(— ) , C 2m _i(— ) = «m + C 2m _ 2 (— ) . 

Then it is easy to verify that for l^Pnf^ ® | — an d v^o = a o = 

|C„(-)| < 2 max(|«!|, |« 2 |) =: M , (8) 

yielding 4 5) (n|-) = sinC„(-), cr^ 5) (n|-) = -cosC n (-). 

From any initial state, \ipo) — a |</?o ) ® 1+/"^ + a< '~' > lv 9 o ) ® l — ) > we then obtain 
at step n 
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and, observing the orthogonality of the |±)^, 

-W(n) = |a«|Vf(n|+) + |a«| 2 <xf \n\-) . (9) 



(T, 



This trajectory of the Turing-head S represents a non-orthogonal pure-state decomposition. 
By using fl7|), (|), (g) (with fl W = = l/y/2) we finally get for \i(t ) = | - Vf& ® | - 1)« 
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■"F^ 77 -)) ^i^l™)) = cos/3 n - (sin^4 n , — cos^4 n ) , (10) 



where (C n (+) - M)/2 < A n = (C 2m (+) + C 2m (-))/2, B n = (C 2m (+) - C 2m (-))/2 < 

(C 2m (+) + M)/2 , n = 2m or 2m — 1. Thus the expression ([K]) indicates that for the local 

dynamics of the Turing head in the 'non-classical' regime the cumulative control loss due 

of of 

to any small perturbation 5 of the given ct^ , oci grows at most linearly with n, so that all 
periodic orbits on the plane jo, <J^\ are stable (see figure ||a - c), as in the case of the 
'regular' control a m = a>\ [0]. This may be contrasted with the chaotic Fibonacci-rule (cf), 
= a m + a m-i (Lyapunov exponent: ln(l + y/E)/2 > 0) , which can be interpreted 
as a temporal random (chaotic) analogue to one-dimensional chaotic potentials ||: each 
step a m is controlled by the cumulative information of the two previous steps. For a small 
perturbation of the initial phase angle «o the cumulative angles A m , B m , respectively, grow 

exponentially with m, and so do the deviation terms AC 2 £(±) = C 2 £'(±) — C 2n p°(±) from 
the periodic orbits (po). Thus the total cumulative control loss induced by the perturbation 
can show chaotic quantum behaviour on the Turing head 0. 

For the Thue- Morse control (|4]) we easily find that for \tpo) = | — Vf s ^ <g> | — lf 1 ^ and 
n = 8m 

C n {+) = 2(a 1 + a 2 )m, C n {-) = 0, 

respectively, which is very similar to the result of the 'regular' machine with a m = (ai + 
res 

a 2 )/2, implying C 8m (+) = 2{ot.\ + a 2 ) and C 8m (— ) = 0. In all cases considered we thus 
find characteristic local invariants with respect to the Turing head (figure H<f). 



3 Parameter sensitivity 

The distance between density operators, p and p', defined by the so-called Bures metric 


D 2 PP > :=Tr{(p-p') 2 }- 

lies, independent of the dimension of the Liouville space, between and 2. For pure states 
we can rewrite 

D 2 = 2(1-|<V#'> T) = 2(1-0') 

of : = KMWfrw umum 2 , 

where the perturbed unitary evolution, U(S), connects the initial state, \ipo(8)), and \ip'). 
This metric can be applied likewise to the total-network-state space or any subspace. In 
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any case it is a convenient additional means to characterize various QTMs: for the regular 

1*6 £C 

case, a m = a% (Lyapunov exponent = 0), and any initial perturbation 5 for p the dis- 
tance remains almost constant ||; for the chaotic Fibonacci rule (cf), on the other hand, 
(®m(p) = «m 5 a m(p') = oi m (p) + 5^ , where 8^ is the cumulative perturbation of the an- 

gle oq^ at step n = 2m — 1) we obtain an initial exponential sensitivity 0. In the case of 
the present substitution sequences we observe for a small perturbation of the given ai,a^ 
no initial exponential sensitivity in the evolution of D 2 , which confirms that any periodic 
orbit is stable (figure |3]a). Finally we display the evolution of D 2 for the total network state 
\if) n ), which also shows no exponential sensitivity (figure [3^, cf. figure 3c in [[|). The re- 
spective distances for tape-spin 1 are similar to those shown. The corresponding behaviour 
under the Thue-Morse control is qualitatively the same. This parameter-sensitivity ]TT| 
has been proposed as a measure to distinguish quantum chaos from regular quantum dy- 
namics. From the results of the present analysis and those in |9| satisfying this criterion we 
conclude that, indeed, only classical chaotic input makes the quantum dynamics in QTM 
architectures chaotic, too (figure H). 



4 Summary 

In conclusion, we have studied the quantum dynamics of a small QTM driven by substitu- 
tion sequences based on a decoherence-free Hamiltonian. As quantum features we utilized 
the superposition principle and the physics of entanglement. Quantum dynamics mani- 
fests itself in the superposition and entanglement of a pair of 'classical' (i.e. disentangled) 
state-sequences. The generalized QTM under this kind of control connects two fields of 
much recent interest, quantum computation and motion in one-dimensional structures with 
'deterministically aperiodic' potential distributions. No chaotic quantum dynamics results 
in this CclSG, clS shown by the lack of exponential parameter-sensitivity. Local invariants 
leading to one-dimensional point manifolds (patterns) exist for a± = ai only. For ol\ ^ ai a 
continuous destruction of these patterns sets in (figure |2|b; hardly visible yet in figure |2|a). 
This reminds us of the disappearing of KAM tori in the classical phase space resulting 
from a small perturbation (see e.g. |12|). Patterns in reduced Bloch-planes jo, , 03 j 
(a quantum version of a Poincare-cut) should thus be similarly useful to characterize quan- 
tum dynamics in a broad class of quantum networks. Due to the entanglement, we can 
see regular, chaotic, and intermediate quantum dynamics, respectively. Furthermore, the 
parameter sensitivity gives a sensitive criterion for testing induced quantum chaos in a pure 
quantum regime. This might be contrasted with the usual quantum chaology, which is con- 
cerned essentially only with semiclassical spectrum analysis of classically chaotic systems 
(e.g. level spacing, spectral rigidity) (T^j. It is expected that a QTM architecture with 
a larger number of pseudo-spins on the Turing tape would still exhibit the same type of 
dynamical behaviour under the corresponding driving conditions. 
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Figure ^ Input-output-scheme of our quantum Turing machine (QTM). 

Figure |]: Turing-head patterns {0,a 2 (n),a 3 (n)} for initial state \ip ) = \ — l/ 5 )® 
| — l^ 1 ) under the control of substitution sequences: (a) quasi-periodic Fibonacci (qf) 
with ol\ = §7T , a 2 = ai + 0.00057T , (6) as in a) but for a 2 = ol\ + 0.037T , (c) as in a) but 
for a 2 — oli + 0.057T ; (d) Thue-Morse (tm) control with a± = |7r , ct 2 = «i + 0.10017T . For 
each simulation the total step number is n = 10000. 

Figure [| Evolution of the (squared) distance D 2 ppl between the perturbed, p' , and the 
reference QTM state, p , under the quasi-periodic Fibonacci (qf ) control: (a) for Turing 
head; (6) for total network state \ip n ) . For p we take \ipo) = \ — l)( s > <8> | — l/ 1 ^ and 
Ql = |vr , a 2 = ocx + 0.03vr . Line A: \%) = exp (-ia[ S) 8/2) |V>o) for p' , 5 = 0.001. Line B: 
\tf Q ) = |^ ) , but a[ = at + O.OOItt , a' 2 = a 2 + O.OOItt for (p') . 
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